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Cn Abstract. The phase space for a system of n qubits is a discrete giid of 2" x 2" points, 

whose axes are labeled in terms of the elements of the finite field F2n to endow it with proper 
geometrical properties. We analyze the representation of graph states in that phase space, 
showing that these states can be identified with a class of nonsingular curves. We provide an 
algebraic representation of the most relevant quantum operations acting on these states and 
discuss the advantages of this approach. 
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1. Introduction 



^ Graph states represent a versatile class of entangled states of uttermost importance in 

quantum information. Originally proposed for implementing measurement-based quantum 

computation [3] ID |5l, they have found numerous applications in other problems 

04 such as quantum error correction lIS^I, secure quantum communication [|7J, entanglement 

^ purification HI, and to demonstrate fractional braiding statistics of anyons ^. Furthermore, 

instances of this family, such as Greenberger-Horne-Zeilinger (GHZ) states and cluster states, 

play an essential role in fundamental tests of quantum nonlocality fTO'TT). Consequently, a 

great deal of effort has been devoted to theoretically understand their properties 1 12 , 13 1 and 

1^ to create and manipulate them experimentally lfT4l [T5l [161 [TTl . 

("^ In addition to all these applications, the graph-state formalism is a useful abstraction 

f^ that permits, in principle, a detailed (although not exhaustive) classification of n-qubit 

^-H entanglement ifTSl [T9l l20l . However, the number of graphs grows very fast with n, so such a 

'^ job becomes rather involved for many qubits. Moreover, for these high-dimensional systems, 

the action of nonlocal operations is difficult to interpret. 

An important remark in this respect is that other quantum tasks (such as, e.g., 
i^ tomography) can be most efficiently implemented in phase space. For n qubits, the phase 

space is a discrete grid of 2" x 2" points 11211 I22ll23ll24ll25]| . whose axes are labelled by the 
discrete Galois field F?" in order to preserve some intuitive geometrical properties [26 1. In this 
way, a bona fide Wigner function can be introduced l! 27ll28ll29l[30ir31J . One natural question 
that comes to mind is whether this phase-space picture may be of any help in reinterpreting 
complex issues of graph states. Our work provides a positive answer: graph states turn out 
to be nothing but discrete curves in phase space 1,32] . with properties closely related to other 
intriguing notions such as unbiasedness ll33l . In this way, quantum operations appear as 
transformations on these curves, with very simple and elegant properties. 
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2. Curves in phase space 

A qubit is realized as a state in a two-dimensional Hilbert space. It is customary to choose 
two normalized orthogonal states, {|0), |1)}, to serve as a computational basis. The unitary 
operators 

CT, = |0)(0|-|1)(1|, CT,= |0)(1| + |1)(0|, (2.1) 

generate the Pauli group ^i of a single qubit under matrix multiplication 1*341. The elements 
of this group are known as Pauli operators and provide a basis of unitary operators on the 
Hilbert space. 

For n qubits, a compact way of labeling both states and elements of the corresponding 
Pauli group ^„ consists in using the finite field F2« (the reader interested in more 
mathematical details is referred, to the excellent monograph by Lidl and Niederreiter ll26l ). 
This can be considered as a linear space spanned by an abstract basis {0i , . . . , 9„}, so that 
given a field element a (henceforth, field elements will be denoted by Greek letters) the 
expansion 

n 

a = ^ fl/ 0; , fl; e Z2 , (2.2) 

1=1 

allows us the identification a <^ (fli,...,fl„). Moreover, the basis can be chosen to be 
orthonormal with respect to the trace operation (the self-dual basis); that is, 

tr(ft0;) = 5,;, (2.3) 

where tr(a) — a + a^ + . . . + a^ , which actually maps ¥2" H^ Z2. In this way, we associate 
each qubit with a particular element of the self-dual basis: qubit , <^ 0,. 

We denote by \a), with a S Fys an orthonormal basis in the Hilbert space of the system. 
Operationally, the elements of the basis can be labeled by powers of a primitive element 
(i.e., a root of a minimal irreducible polynomial, called the primitive polynomial), and reads 
{|0), |(7), . . . , |c7^ ^' = 1)}. These vectors are eigenvectors of the operators Zp belonging to 
the generalized Pauli group 3^„, whose generators are 

a a 

Here the additive characters x are defined as x(a) = exp[/7rtr(a)]. We have then 

ZaXp=X{al5)XpZa, (2.5) 

which is the discrete counterpart of the Heisenberg-Weyl algebra for continuous variables. 



The operators ( 2.4 1 can be factorized into tensor products of powers of single-particle 
Pauli operators. This factorization can be carried out by mapping each element of ¥2" onto an 
ordered set of natural numbers according to 

Z„ = afi «) • • • (g) o-f , Xp = of' (g) • • • ® cxf" , (2.6) 

where a,- = tr(a0,) and bi — tr(j3 0,) are the corresponding expansion coefficients for a and j3 
in the self-dual basis. 

The commutator of two monomials in ^„ is 

[Za.Xp^ ,Za,Xp^] = [Z(aij32 + jSiaz) - l]Za,Xp^Za,Xp^ , (2.7) 

so they commute when 

tr(aii32 + i3ia2)=0. (2.8) 
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We next recall ETl l22l |231 that the grid defining the phase space for n qubits can be 
appropriately labeled by the discrete points (a,/3), which are precisely the indices of the 
operators Za and Xn: a is the "horizontal" axis and )3 the "vertical" one. 

From this perspective it seems reasonable to consider a (parametric) curve in this discrete 
grid as a set of 2" + 1 points (a(K'),)3 (fc)), where the parameter K runs through the field F2'>. 
The curve pass through the origin when (a(0), j3(0)) = (0,0), and it is nonsingular (i.e., with 
no self- intersection) if all the pairs (a(K'),j3(K')) are different. 

If we represent the curve as 

a{K)J'fa.K^\ P{K) = Y^PiK^\ (2.9) 

/=() 1=0 

with a,-,j3; e F2", we will say that the curve is additive commutative whenever 

£tr(a,i3;)=0. (2.10) 

Due to the form in equation ( |2.9| l, the following property is automatically fulfilled 

a{K + K') = a{K) + aiK'), I3{k + k') ^ I3{k) + I3{k'), (2.11) 

which means that by summing the coordinates of any two points of the curve we obtain 
another point on the curve. The condition ( |2.10[ ) guarantees that the monomials labeled with 
points of such a curve form indeed an Abelian subgroup under multiplication [notice the 
condition p.8| l]; i. e., 

[Za{K)Xp{K),Za(K')Xli{K')] =0. (2.12) 

(Br 
2 



Observe also (2.11 1 implies that this group is isomorphous to Z®" 



The 2" monomials {Za(K)XRiK)} constitute thus the whole set of stabilizer operators ll34l 
and n of them appropriately chosen can be considered as generators. This means that one can 
pick n specific points that generate the whole curve by simple addition. The simplest way 
is to pick the horizontal coordinates as the points corresponding to the elements of the basis 
a ~ {9i,...,9„}, the order of points being unessential. 

The simplest form of additive commutative curves are the straight lines 

a{K)^IXK, j3(K-) = VK-, (2.13) 

which can be represented in the regular form a = 0,j3 = Aa. It is a well established 

result fSS", '361 that the operators {ZaXa^xa] associated to straight lines commute for any 

fixed value of A G F2", while the eigenstates of the set {Za} define the standard logical basis. 

A curve is called nondegenerate (or regular) when can be represented in the explicit form 

i3 =/(«) = '£'/''«"' or a = ^(j3) = X;v^,i32', (2.14) 



with 0;, I//,- e F2n . Such curves are nonsingular and the commutativity condition ( 2. 10 1 imposes 
the restrictions 

<Pk = <P^-k. Wk^wtk, ^=l,...,[(«-l)/2], (2.15) 

where [] denotes the integer part. For even values of «, the additional requirement (^„ « = (^^/j 

9"/- 

iWn/2 — W„/2 ) should be satisfied. 

The expansion coefficients are related to the curve equation in a very simple form. 
Indeed, by rewriting the equation j3 — f{a) as 

j3=/(a)=/(f tr(a0,)0,Uf tr(a0,)/(0,)-f X;'a2' (2.I6) 

\k=l ) k=\ k=\ i=0 
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and comparing with ( |2.14| i, we immediately get that 

0. = LC/(0O- (2.17) 

The degenerate curves have more involved structure ll32l . The degeneration basically 
means that a and j3 do not take some values in the field, which is compensated by multiple 
appearance of other (admissible) "coordinates" , since the curve is nonsingular Such curves 
are characterized by the degrees of degeneration r^ and rj^ . This degree means the following: 
if {aj,l5j) is a point of a degenerate curve, for each aj there are 2"^'P values of )3, such that 
the points {aj.fik) (^ = I7 • ■ ■ ,2"^'^P) belong to the same curve and, conversely, for each j3y 
there are 2"^'^° values of a, such that the points {ak,Pj) also belong to the same curve. 

3. Operations with curves 

We next show how the action of Clifford operations can be described in this finite grid. For a 
single qubit, the simplest of such operations are x- and z-rotations, which are represented by 
the unitaries Ux^z = exp(i;r(7r ,;/4). They produce the transformations transformations 

X— rotation : {c7^ 1-^ Gx, Oy i-^- a,, a, i-^- Oy] , 

(3.1) 

z— rotation : {(7^ i-^- a,., 0|y i^ <7v, c- 1-> cTj.} . 

Multiqubit rotations can be labeled by ^ e F2" , so that in the self-dual basis the expansion 
^ = Y^i £,i9i indicates on which qubit the rotation is performed: if ^,- = 1, then the i-th qubit is 
rotated. The general form of a local x- and z-rotation is thus Ux,z = 11"= 1 "x.z ■ The indices 
of the monomials ZaXg transform then according to 



X— rotation : 



z— rotation : 



i3^j3+lA.^itr(a0,)0,, 

(3.2) 



With local transformations it is always possible to transform any degenerate curve into a 
nondegenerate one. Moreover, any curve can be cast into the form /3 = /(a), where f{a) is 
an invertible function, so that the curve equation can be also written as a = /^' (j3). 

One- and two-qubit gates are crucial in quantum computation. As an archetypal example, 
we consider the nonlocal XOR,j operator (applied to the /th and yth qubits) with the following 
action 

1 ('■) di^'' ^ a^ di''' , ai'^ 1 f^') ^ ai 1 <'' , 

(3.3) 

In our algebraic language this can be elegantly represented by 

aH^a + 0,tr(a0;), j3 1-^ j3 + 0^tr(j30,), (3.4) 

so that 

xOR,;=^|A+tr(A0,)0i)(A|. (3.5) 
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Table 1. IiTeducible polynomials and self-dual bases for the Galois fields F2" used in this 
paper. 



Field Irreducible polynomial Self-dual basis 



% 


x2+.v-+l = 


{a,G^} 


F23 


a:3+a:+1 = 


{a\a\(j^} 


F24 


x*+x+\=0 


{CT^CT^CT'2 cr"} 


F25 


x'+jc^ + l =0 


{CT3,CT',CT",fTl5 ff"} 



In particular, the curve j3 = f{a) is transformed into 

/3 ^ /(a) +tr(a0,)/(0,) +tr[/(a)0,] 0; +tr(a0,) tr[/(0,)0/]0, ■ (3.6) 

Equivalently, the expansion coefficients are ^k '-^ ^k + &} f{^i) +/^ (^z) ^j ■ For example, 
for two qubits, the curve )3 = a is transformed, under the action XOR12, into j3 = a^a; where 
(7 is a primitive element of F22 , solution of the irreducible polynomial x^ + + 1 =0. For three 
qubits, the curve j3 = a becomes, also under the action XORn, j3 = aa + a^ + a^, where 
<7 is the primitive element for Fjs, solution of jc^ +x+ 1 = 0. In Table n] we summarize the 
relevant information for the Galois fields used in this paper. 

The SWAP operator exchanges the states of the ith and 7th qubits; i.e., 

SWAP;^|...,fl,-,...,a^-,...) = |...,aj,...,a,-,...). (3.7) 

At the level of field elements, its action is 

aH^a + etr(ae), j3 1-^ j3 + etr(/3e), (3.8) 

where e = 0, + Oj. Consequently, we can write 

swAP,7 = ^|A+etr(£A))(A|, (3.9) 

A 

so it transforms the curve j3 = f{a) into 

j3 ^ f{a) +tr{ae)f{e) + tr[f{a)e] £ + tr(ae) ti-[/(£)£] e . (3.10) 

The expansion coefficients of the transformed curve are (j)k + £^ /(£) + /^ (£) £. 
As our final nonlocal operation, we take squeezing, represented by 

S^=J^\^X}{X\, (3.11) 

A 

whose action on the basic monomials is S^ZaXaSV^ ~ Zf„Xc_io, so that, the curve is 

transformed into j3 = ^^^f{a), which, as expected, is just a global scaling 

One can also define curve composition in such a way that given two curves j3 = /(a) 

andj3 =g(a), the function j3 =f{g{a)) is also a curve iff /(g(a)) ^g{f{a)). 

The last remark concerning these curves is that they can always transformed into the 

horizontal axis j3 = 0. Indeed, one can always find an operator Pf such that 

PfZaPf' - (p/(a)Z„X^(„) , (3.12) 

where (Pf{a) is a phase factor The operator Pf can be conveniently expanded in the the basis 
{|A)} of eigenstates of X;^ 

^/ = E4^'|A)(A|. (3.13) 
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If) 
The coefficientsc^ can be determined by noting that 

PfZaPj' =l^cf_^J['*\Xra){l\ 



(3.14) 



while 



Za^/-(«)=Lx[A/(a)]|A + a>(A|, 



(3.15) 



Sf) 



so that cl satisfy the recurrence relation 



Af) If)* 



,Pf{a)x[Xf{a)]. 



(3.16) 

The phase ^f{oc) is evaluated by substituting A = o: in the above equation and taking Cg = 1, 
which finally leads to 



'^T4^=xWfi^)]c[%. 



(3.17) 



„(/) 



The solution of equation (3.17 1 is not unique, but choosing all Cg for a basis to be positive, 
we arrive at the following explicit expression 



Jf) 



nv(-i)^'^^ 



(-1) 



£'M£ 



(3.18) 



Here £' = 



.,£„) is the vector of the expansion coefficients of |A) in the self-dual basis. 



(/) 



and the principal branch of the square root is chosen. In addition, Mfj = tr[a,/((J/))] and 



M: 



M 







(/) 



J <h 



j>i. 



(3.19) 



All this means that given a curve we can immediately obtain the eigenstates |A;/) of the set 
of commuting monomials corresponding to this curve: |A;/) — f/|A), with |A) being the 
computational basis. 

In a similar way, the curves a = g(j3) can be transformed into the vertical axis a — 0: 

Zg(p)Xp^xifiP)QgXpQg\ (3.20) 

where y/{P) is a phase. Again, one can show that Qg can be expressed as 

e, = £cif'|J)(X|, (3.21) 



(g) 



and the coefficients c)^ satisfy (3.17 ). 



4. Graph states and their phase space representation 

Consider a graph G — {V,E}, which is a set of vertices V = {1 , . . . , «} connected in a specific 
way by a set of edges E cV xV, specifying the neighborhood relation between vertices 1 37 1 . 
In the following, we only consider simple graphs, containing neither loops neither multiple 
edges. To each vertex ; we attach a single qubit and each edge {/, j} represents the interaction 
between the corresponding qubits. With the graph we associate a complete set of n commuting 
monomials 






(4.1) 
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where A^, is the neighborhood for that vertex, A^, = {/' e y|{/,y} e E} (i.e., all the vertices 
connected to /). Please, note that we have interchanged the order of (7^ and Ox with respect to 
the standard definition fT2l|: this corresponds to a global Clifford y-rotation, and will allow us 
to work with curves j3 = f{cc) and not a — g(j3). 

The associated graph state \G) is just the unique common eigenvector of {Ki} with all 
eigenvalues +1. The monomials {Ki} are also the generators of the stabilizer of \G) 1*381, 
which is an Abelian subgroup of the local Pauli group ^„ |34|. Any stabilizer state is locally 
equivalent to some graph state ||39l (more exactly, they are equivalent under local Clifford 
transformations). 

The crucial observation for what follows is that the graph-state generators {K,} can be 
written precisely as {ZaXn^^^^'j} , when a runs the self-dual basis {9i, ..., On}, if the function 
/3 = /(o:) satisfies the extra condition 

ti-[0,/(0,)]=O. (4.2) 

Note in passing that the .x-rotation U^ , with 

<^=^tr[/(0,)0,]0M (4.3) 

reduces any curve /3 = f{cc) to the form /3 = /(ot) +£,tr[/(0,)0,] tr(a0,) 0;, which obviously 



satisfies the graph-state condition (4.2 1. For instance, in F23, the non-graph curve j3 — a a is 
transformed into the graph curve j3 = (J^a^ + CJ^a^ by two x-rotations in qubits corresponding 



■J c 3 5 

to elements of the self-dual basis a and a . u" m? . 



It is indeed possible to find a generic form of the curves that fulfill (4.2 1. To this end, let 



us substitute the explicit expansion (2.14 1 into the condition (4.2 1; 

-1 



^tr(</.,-02'0,)=O, k^l,...,n. (4.4) 



Since the expansion coefficients 0, can be expressed in the form of equation (2.17 1, the 
condition tr(0o) — should be always satisfied for graph-state curves. When « is odd, and by 
using the property ( |2.15| l, we can reduce the sum ( |4.4| l to the simpler form 

tr((^2'"'0,)=O, (4.5) 

which shows that any curve for an odd number of qubits with ^0 = corresponds to a graph. 
When n is even, we have to take into account an additional condition on the coefficients, 
^n/2 = 0^/2 ' which leads to the restriction 



(/)o = £tr((/.„/20r^'0,)02. (4.6) 



(=1 

A graph can be also characterized by its adjacency matrix F, which is a symmetric nxn 
matrix such that 

r,/ = <^ (4.7) 

I 0, {i,j}(^E. 

The adjacency matrix for graphs associated to the curves fulfilling ( |4.2[ ) has the form 

F„=tr[0,/(0,)]. (4.8) 



For an arbitrary curve j3 = /(«), the result still holds once the diagonal part of F,y in (4.8 1 is 



removed. Note, however, that the same adjacency matrix may correspond to different curves. 
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For instance, for three qubits the curves on F23 j3 = (7'*a and j3 = (7'*a^ + iJ^a'^ have the 
same adjacency matrix: the second curve corresponds to a graph, while the first one does not. 
Given the matrix F of a graph G, one can use (4.8 1 to immediately find that the 
corresponding f{9j) are 



/(0,)-£ti-[0,/(0,)]0, = £F,,-0,- 



(4.9) 



!=1 



Multiplying ( |4.9[ ) by Uj = tr(a0y) and summing we obtain 

f ajfiOj) = f fiajOj) - /(a) = f a;F,-,0,. , 
7=1 )=i 'j=i 

so the curve corresponding to the matrix F is 

j3 = f tr(a0,)F,-,0,- = £ E To^f 0,«2\ 

ij=l k=li,j=l 

or, equivalently, its expansion coefficients 0^ are 

u=i 

Let us work out some examples. The simplest one is perhaps a linear chain of n qubits, 
with adjacency matrix 

/OlO 00\ 

10 1 ■■■ ■■■ 



(4.10) 



(4.11) 



(4.12) 



10 1 
: ■■ ■■ 



■■ 
■■ : 



(4.13) 



0:0 1 10 
VOOO 1 00/ 

so the associated curve is 

n-l 

/3 = E tr(a0,)(0/+i + 0,-1) +tr(a0i) 02 + tr(a0„) 0„_i . 



(4.14) 



Another interesting example is a star graph, in which every single vertex is connected only to 
the center vertex (which we take as the nth qubit). Now the adjacency matrix is 

/ ... 1 \ 









Vl 1 
which leads to the very simple curve 

j3 = 0„ + tr(a0„). 



1 

1 

1 0/ 



(4.15) 



(4.16) 
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-f r 



n-l 



/? = ^tr(a^i)(^i+i+^.- 



i-i 



i=2 



+ tr(afli) O2 + tT{aQn) 0„ 





= 6,, + trM„) /? = 5^ ti{a6i+i){6,+2 + ^i) 

4=1 

Figure 1. (Color online) Curves associated with different graph states. In all the cases, we 
have restricted ourselves to five qubits, although the curve equations are written for an arbitrary 
number n. 



Finally, we consider a graph consisting of a closed polygon, with matrix 
/ 1 ••• 1 \ 

10 1 ■■■ 

1 ■. ■. ■■■ : 
: ■. ■■. 10 




Vl 













1 



1 

0/ 



(4.17) 



and associated curve 



j3 = ^tr(a0,+i)(0,+2 + 0O, 
(=1 



(4.18) 



where the indices of the basis elements are taken mod n. 

Please, observe carefully that for each of the examples worked out so far the following 
relation holds: 

/3(0,)=E^i' (4.19) 

that is, the curve evaluated at each vertex appears as a sum of basis elements associated with 
all the connected vertices. In figure 1 we graphically summarize these relevant examples of 
the correspondence between phase-space curves and graphs. 

Concerning the nonlocal operations XOR and SWAP, one can check that graph curves are 
transformed into graph curves. 

Local Clifford operations on graph states can be easily explained in terms of the local 
complementation. For the graph G = {VjE}, the local complement of G at / is obtained by 
complementing the subgraph of G induced by the neighborhood A^, of / and leaving the rest of 
the graph unchanged. This is implemented by the local unitary 

V, ^4U 4 (4.20) 
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The curve coefficients transform thus in a simple way 

h ^ h+f+'m + 1 ti-[/(eO0<] 0f+' . (4.21) 

' (=1 

It is clear form that the local complement does not change the coefficient 0o, which reflects 

the fact that a graph is also transformed into a graph. 

(f) 
Finally, it is worth noting that for graph curves, the coefficients cV in the expansion 



(3. 17 1 simplify to 



c^P = i-l)^''" , (4.22) 

where F is the adjacency matrix corresponding to the curve. 

5. Factorization structure 

One essential property of these curves is their factorization structure, which is related to the 
commutation condition between blocks of single-qubit operators. 

Consider the commuting monomials labeled by points of a given curve. Each monomial 
is a direct product of n Pauli operators. Let us divide each monomial into two parts, so that 
the first part contains k Pauli operators, corresponding to the (/i , . . . , 4) qubits and the second 
part n — k operators of the rest of the qubits. If any Pauli operator in the first block commutes 
with all the others in that block, we will say that the corresponding set of monomials is 
factorized at least into two subsets. Obviously, the second blocks would then also commute 
among themselves. Moreover, inside the first or second blocks some sub-blocks may exist 
that commute with corresponding sub-blocks, etc. Thus, we can represent any curve 7 e ¥2" 
in the following factorized form: 

7=={mi,m2,...,mAr}, (5.1) 

where < mi < m2 < . . . < niN and m, e N is the number of particles in the /th block that 
cannot be factorized anymore. It is clear that {mi , m2, . . . , »%} is just a partition of the integer 
n, so the maximum number of terms is n, which corresponds to a completely factorized curve, 
7= {1,1,. ..,1}, and the minimum number of terms is one, corresponding to a completely 

nonfactorized curve 7 = {«}. 

The rays (straight lines passing through the origin) a = 0, j3 = 0, and j3 = a are always 
completely factorized, while all the other are apparently completely nonfactorized for any 
dimension, except the special case of four qubits, where there are two locally equivalent rays 
with factorization (2,2). 

It is easy to convince oneself that completely factorized curves can be written as 

/(a) = X^ti-(a0,)^,-0,, (5.2) 

1=1 

where <^, are the components of an arbitrary £, £ F2" in the self-dual basis. 

For bifactorized curves, whose adjacency matrix is just a 2k x 2k permutation matrix, the 
corresponding form is 

f{a) = Y^ [tr(a0,„ ) 0,„ + tr(a0,„ ) 0,„] , (5.3) 

where {ii,ji},{i2,J2},---,{^kijk} are ordered pairs of connected vertices. In this case, two 
relations are automatically satisfied; f[f[a)) — a and tr[/(a)] — tr(a). 
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{1,1,1} 

{2,1} 

{3} 



O P=o''a 



X_J |3=aa 



O O O P=o'^a 




Figure 2. (Color online) Inequivalent curves for three (top) and four (bottom) qubits. In the 
left, we show the factorization structure; in the right, the curve in the simplest algebraic form 
and in the middle the locally equivalent graphs. 



In general the problem of fully determining the factorization of a curve is quite difficult, 
and basically corresponds to the separation of an adjacency matrix in independent blocks by 
permutation of lines and columns. The number of such blocks can be determined by, e.g., 
counting the number of zero eigenvalues of the Laplacian matrix |40|. Nevertheless, some 
general considerations can be made. Let us evaluate a function /(oc), corresponding to a 
graph, in an arbitrary element 0/ of a self-dual basis. 



m) 



:^tr[/(0,)0;]0; 






Qj: 



(5.4) 



which is the sum of elements of the basis corresponding to vertices (qubits) connected to 
the /th vertex. This means that this /th vertex is disconnected form the rest of the graph if 
f{Oi) =0. Moreover, for each disconnected subgraph A,, the function /(a) maps the subgraph 
onto itself /(A,) —Aj and thus 



tr/ E0; 







(5.5) 



\j(^i 



for each separated block. 

In figure 2 we list the locally inequivalent curves for three and four qubits, where they 
are expressed in the simplest algebraic form (not necessarily corresponding to graphs) and the 
equivalent graphs are plotted. This can be easily continued to higher number of qubits. 



6. Some applications 

6.1. Local and nonlocal transformations 

The algebraic properties of the commuting stabilizers {ZaXf/(^\} allows us to scrutinize in 
a very simple way the action of both local and nonlocal transformations. Moreover, given 
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then relation between curves and adjacency matrices we can easily find how graphs are 
transformed, avoiding in this way involved procedures |39|. 

Let us start by considering the linear operations treated in section 3. It is clear from 



(3.3 I that x-rotations (labelled by the parameter t^, whose nonzero components indicate the 
qubits subjected to rotations) do not change the graph; they add only a diagonal element to 
the adj acency matrix: F i-^ F + diag (i^ i , . . . , ^„ ) . 

On the other hand, z-rotations change the adjacency matrix in a nontrivial way, as 



hinted by (3.3i. First, we find the corresponding curve for which we change the variable 
a t-^ jj. = (X + Y,"=i ^(tr[/(a)0,] 0,-. This requires the inversion of the system 



^, = a, + ^,-^ayF, 



U' 



(6.1) 



.7=1 



(6.2) 



to be able to determine a = a(/i) and thus find the set {Z^X„(„\}, where g = /(o:(/i)). 
Finally, the curve j3 = g{lJ.) can be converted into a graph-curve (if necessary) by some x- 
rotations, which immediately gives the new adjacency matrix. 

As an example, let us analyze a four-qubit graph with adjacency matrix 

/ 1 ^ 

and find how it is transformed un der z -rotations of the second and the fourth qubits. The 
curve corresponding to the matrix (6.2i is j3 = a^a^ + (j''a^. The rotation 1 (g) m, (g) 1 (g) m, 
transforms it into the curve j3 = (7 ' ^a + a 

1 





1 







1 



(jl5«4_ 






This means that 



F' 



1 




1 

1 J 



(6.3) 



The final adjacency matrix is obtained from the above matrix just by removing the diagonal 
elements (which can be done just by applying the jc-rotation 1 (g) 1 (g) 1 ® Ux). The resulting 
graph-curve is j3 = a^a^ + a^'^a'^ + a^a^. 

Next, we turn our attention to nonlocal operations. According to (|3.6|l the adjacency 



matrix associated with the graph-curve transformed by XOR,j is 

-tr[0,/(0,)]5//+tr[0,/(0,)]5; 



'kj- 



(6.4) 



F„-^F,,=tr[0,/(0,)]- 

This can be rewritten as F i-^- F' = GTG^, where Gu = 5ki + Skjdu is a triangular matrix 
with det G =\. For instance, the application of the gate XOR12 to the graph defined by the 
matrix ( |6.2| ) leads to 

1 




F' = 



1 
1 

VO 



(6.5) 



0/ 



C7^a^. 



and the corresponding curve is )3 

Finally, the squeezing operation ( |3.11| l transforms the curve j3 = /(o:) into j3 
/(/i^^')^^' and, thus the adjacency matrix F corresponding to the initial graph into 



-,j ^ tr[0,-/(0,-^ 



-1^ 



ri=tr 



Ltr(0,^ 

k=l 



-1 



ek]ekf 



Ltr(0,-^ 

w=i 



-1 



B,)B, 
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= £ ti-(0,^-'0,)tr[0,/(0,)]ti-(0,-^-'0/)=M,,r„M/,-, (6.6) 

k,l=\ 

where M,vt = tr(0,^ " ' %) is the mattix M for the ray j3 = ^ " ' a. 

For example, the graph corresponding to the adjacency matrix (6.2 1 is transformed under 
Sa (where a is the primitive element for F24) converts into the graph defined by 

/ 1 1 1 \ 
r' 10 1 

r = 10 ' ^^-^^ 

\ 1 1 0/ 
and the corresponding curve is j3 = a'^a^ + a^a^. 

6.2. Computing reduced density matrices 

As final illustration of the potential of our method, we consider the density matrix of the graph 
state corresponding to the curve j3 = /(o:), which is given by 

p=2-"£z„X^(„). (6.8) 

a 

Now assume that we trace over n — k qubits. This leads to the following reduced density 
matrix 

P=2" Y, Zaiei+...+atetXf{a^ei+...+aM ■ (6-9) 

With the machinery developed so far, this can be recast in the explicit form 

p^2-' £ j^^^a,^tr[/(a,a,+...+a,a,)a,]^2-*- £ n^af'aF'^' ""^'" , (6.10) 
«! aj.GZ2 i=l ai....,af.£Z2 i=l 

where F is the adjacency matrix of the initial graph. 

7. Conclusions 

In summary, we have developed a phase-space picture of stabilizer states, i.e., a complete 
set of commuting monomials in the generalized Pauli group. This alternative description 
identifies these states with functions (which can be associated with specific nonsingular curves 
in discrete phase-space) over the finite field Fj. 

All the operations on graph states can be expressed in this language as transformations 
of the curves. Some of these transformations, such as the XOR gate or the discrete squeezing, 
are simpler to apply to the algebraic functions than to the graphs. On the contrary, the 
factorization is more involved at the algebraic level. We expect that the proposed method 
can be used for classification of graph states when the number of qubits is large and graphs 
do not provide any useful information. Work in this direction is in progress. 
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